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Abstract: We report on the observation of single colour centers in 
natural diamond samples emitting in the near infrared region when optically 
excited. Photoluminescence of these single emitters have several striking 
features, such as a narrow-band fully polarized emission (FWHM 2 nm) 
around 780 nm, a short excited-state lifetime of about 2 ns, and perfect 
photostability at room temperature under our excitation conditions. We 
present a detailed study of their photophysical properties. Development 
of a triggered single-photon source relying on this single colour centre is 
discussed in the prospect of its application to quantum key distribution. 
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1. Introduction 

Significant efforts have been expended recently in the development of new and reliable single- 
photon sources, due to an increasing number of applications of such sources in quantum òptics, 
namely in quantum cryptography and for the realization of quantum gates [ 1 ] . For instance, bà- 
sic security proofs of quantum key distribution (QKD) rely on encoding information on single 
quantum objects, namely single-photon light pulses [2, 3]. Most practical QKD systems use 
weak coherent pulses as an approximation of single-photon sources, which greatly simplifies 
the experimental implementation [4]. However, faint làser pulses with mean photon number 
smaller than unity can still contain more than one photon due to Poissonian photon number 
statistics. Increasing attenuation on the quantum channel to lower the probability of multipho- 
ton pulses is a practical solution, but it also reduces the secret key creation rate. Moreover, due 
to dark count noise of the detectors, the vulnerability of the most commonly used QKD Bennett 
and Brassard BB84 protocol [2] increases with losses on the transmission channel [3]. To im- 
prove security of long-distance QKD, one as either to conceive of new protocols, like the one 
relying on decoy states [5], or to switch to single-photon sources. 

Single photons can be conditionally prepared using a pair emission process like atòmic cas- 
cade [6] or spontaneous parametric down-conversion [7]. One photon of the emitted pair is 
used as a "trigger" which prepares the field modes corresponding to the second photon in a 
single-photon state. With this scheme, heralded single photons at telecom wavelengths can be 
efficiently prepared using bulk^' 2 ' nonlinear crystals [8], or guided-wave nonlinear devices [9]. 
These systems have been experimentally confirmed as promising for optical fiber-based long- 
distance QKD [10] . 



However, emission time remains random for this type of single-photon generation. In quan- 
tum cryptography-based communication system, it is highly desirable to deal with single pho- 
tons synchronized on an external clock. 

Clock-triggered single-photons can be efficiently produced by suited pulsed excitation of 
the emitting dipole [11], which should undergo a full cycle, including excitation, emission and 
reexcitation before it emits the next photon. Since first proposal by De Martini et al. [12] a wide 
variety of schemes to implement a single-photon source have been worked out, relying on the 
fluorescence of an isolated single emitter, such as a molecule, an atom, a colour centre, and a 
semiconductor quantum dot (see Ref.[l] for a review). 

A reliable room-temperature clock-triggered single-photon source was recently realized, 
based on the pulsed optically excited photoluminescence of a single nitrogen-vacancy (NV) 
colour centre in a diamond nanocrystal [13]. This turn-key single-photon source was then im- 
plemented in QKD experiments using BB84 protocol with polarization-encoded single-photon 
sequences [14, 15]. As theoretically expected [3], use of single-photon light pulses showed 
measurable advantages over faint làser pulses in terms of security of the quantum channel and 
key creation rate in the high attenuation regime. 

The wide variety of photostable colour centres in diamond [16] offers high opportunity to 
optimize single-photon properties so that they better meet application requirements for prac- 
tical QKD. For instance, polarized single-photon emission, together with nanosecond-range 
excited state lifetime and nanometer spectral bandwidth, should be achieved. Recently, the 
nickel-related NE8 defect was observed at the single emitter level, both in natural diamond 
sample [17] and in CVD diamond film [18]. This colour centre, which consists in a nickel atom 
in divacancy position and four nitrogen àtoms in its first coordination sphere, is characterized by 
a weak electron-phonon coupling. Most of its photoluminescence is then concentrated within 
its zero-phonon line, corresponding to a sharp emission line (FWHM 1 .5 nm) around 800 nm. 

We report on the observation at room temperature of single diamond colour centres with per- 
fect photostable luminescence at room temperature, in the same wavelength range as the NE8 
colour centre. Assumed to be also nickel-nitrogen related impurities, we study in details their 
photoluminescence properties under optical excitation with a red cw làser diode. Both single- 
photon emission and shelving in a metastable level are observed from photon time-correlation 
measurements, with an Hanbury Brown & Twiss setup. We show that these emitters meet the 
required properties for practical single-photon QKD as stated above. First, a narrow emission 
band is observed around 782 nm, which is both in the open-air and fiber-optics optical commu- 
nication window. Second, an excited-level lifetime around 2 ns is inferred, which is compatible 
with narrow time-window analysis in order to reduce the effect of photodetector dark count 
noise. Third, light emitted is fully linearly polarized, which is perfectly suitable for polariza- 
tion information encoding on the emitted single photons. 

2. Experimental set-up and results 

We use an home-made scanning confocal microscope to select the single defects in the sample 
(see Ref.[19] for a detailed description of the setup). The detection system includes a Hanbury 
Brown & Twiss (HBT) setup with two silicon avalanche photodiodes (APDs) in single-photon 
counting regime (SPCM, Perkin Elmer) on each side of a 50/50 non-polarizing beam splitter. It 
also integrates a spectrograph with a cooled CCD array for recording spectrum. We first raster 
scan the sample (2 x 2 x 0.32 mm type lla natural bulk diamond wafer from Elementó, The 
Netherlands) to identify well isolated photoluminescent emitters using a làser diode emitting at 
687 nm as cw excitation source. The làser beam is focused on the sample about 4 fim below 
its surface by a microscope objective (x 100, NA= 0.95) which is also used for emitted light 
collection. The làser focus spot is about 1/im FWHM. Fig.l(a) displays a 9 x 9 fim scan of the 



sample, showing photoluminescence from a single nickel-nitrogen related defect. The z-axis 
value indicates the count rate of one APD of the HBT setup. From the scan in Fig.l(a), the 
observed signal to background ratio of the photoluminescence is as large as 60 : 1 . 
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Fig. 1. (a) Fluorescence intensity raster scan of a natural diamond sample showing lumi- 
nescence from an isolated colour centre. The counting rate corresponds to the output of one 
detector in the HBT setup. Its maximum value is about 4 x 10 4 counts/s. (b) Photolumines- 
cence spectrum from the single emitter observed in Fig.l(a) with 10 s integration duration. 
The spectrum displayed has been corrected for the quantum efficiency of the CCD which 
varies from 85 to 53% in the 740-840 nm range considered. The narrow peak (1) at 782 nm 
(about 2 nm FWHM) is the zero-phonon line of the nickel-nitrogen related defect [17]. The 
sharp peak (2) at 756 nm is related to the one-phonon Raman scattering line of the diamond 
lattice with 1332 cm -1 frequency shift. Inset: zoom on the zero-phonon line, showing more 
clearly the intensity of the phonon wing. 

The HBT setup is used to investigate the fluorescence intensity correlation function. Thanks 
to the narrow spectral emission, a bandpass filter placed before the HBT setup allows us to se- 
lect the photoluminescence from the defect while rejecting Raman scattered light due to the di- 
amond crystalline lattice. Appropriate data processing of the photodetection events in the HBT 
configuration gives the histogram of the time separations between two consecutively detected 
photons. This delay function is equivalent to the second-order intensity autocorrelation func- 
tion for short time scales [20]. Evidence for single centre emission is given by the observation 
of antibunching [21]. This effect appears as a dip around t = in the recorded delay function 
(see Fig.3(a)), as the consequence of finite-time cycle of excitation-emission-reexcitation for a 
single dipole. 

For each observed single emitter, we simultaneously record the spectrum of the emitted light 
on the spectrograph. Fig.l(b) shows the spectrum of the light collected from Fig.l(a) emitter. 
Narrow emission peak at 782 nm is tentatively associated to a nickel-nitrogen related impurity 
in diamond [16, 22]. From the spectrum, we infer that even at room temperature about 70% of 
the total intensity is concentrated in the zero-phonon line. This value is much higher than that 
of nitrogen-vacancy colour centre, for which it is only of the order of 1%. 

We also investigated the polarization properties of the single colour centre relative to absorp- 
tion and emission of light. We have monitored the photoluminescence intensity while rotating 
the excitation làser linear polarization angle and obtained oscillations with a contrast of about 
96%. This contrast value proves that the single colour centre behaves like a dipole relative to 
absorption of light. Using the quarter wave plate method, we studied the emitted light polar- 
ization properties. We observed that the collected light is perfectly elliptically polarized with 
an aspect ratio of 0.25. This measurement indicates that the single defect also behaves like an 



emitting dipole. The linear polarized light becomes elliptical after propagation though òptics 
including high numèric al aperture objecti ve and dichroic mirror, both known for modifying the 
polarization state of light. 

3. Photophysics of the 782 nm emitting colour centres 

In order to investigate photoluminescence dynamical processes, we model the response of the 
single colour center within the framework of a three-level system corresponding respectively to 
ground level (1), excited level (2), and metastable level (3) as shown in Fig.2. 
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Fig. 2. Three energy levels scheme for the single colour centre emitting at 782 nm central 
wavelength, with corresponding decay and intersystem crossing rates. 




The single colour centre is excited from its ground to excited level. It then either goes back to 
the ground level, or it non-radiatively decays to the metastable level, which presence is asserted 
by photon bunching in the intensity correlation measurement. From this "dark" level, the system 
finally goes back to the ground level by a phosphorescence process. The dynamics can be 
described by a set of rate equations on the populations {pi,p 2 ,P3} of the three considered 
levels 

pi= -r\2P\+r2\p2 + n\p^ 

P2 = r u pi - (r 2 i + r 2 i)p2 (1) 

P3 = r 2 3P2~r3lP3 

where r nm (n,m = 1,2,3) denotes the transition rate from the level («) to level (m). We assume 
that the transition rates r 2 3 and rçi respectively to and from the metastable level are small 
compared to the excited leve to the ground-level decay rate r 2 \. Following Ref. [23], we can 
infer, for the normalized intensity correlation function 

(2)(0 = mm) (2) 

8 V)- {I[t))2 W 
the following analytical expression, vàlid for t > 

g {2 \t) = p 2 (0-t)/ P2 H = 1 - (1 +a)e-^+ae-^, (3) 

Ai = ri2 + r 2 i, (4) 

A2 = ni +r 2 3n 2 /(ri 2 + r 2 i), (5) 

a = n 2 r 2 3/[r3i(rn + r2i)], (6) 



where p 2 {°°) is the steady-state population of the excited level (2), and where p 2 (Q;t) denotes 
the population of level (2) at time t, starting from level (1) at time t = 0. Recalling that r 2 3 and 



r3i are small compared to ri\, it follows that X\ is much larger than Two limit expressions 
of Eq.3 can then be found 

(i) On "short" time scale (t < 20 ns), the g' 2 ' function exhibits photon antibunching: 

g( 2 \t)^\-{\+a)e-^. (7) 

(ii) On "long" time scale (f > 20 ns), the antibunching associated with luminescence of the 
single emitterdoes not affect the g^ function [24] , so that: 

g {2 \t)~\+ae- X *. (8) 
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Fig. 3. (a) Number of photon coincidences c(t) (right scale) recorded with the HBT 
setup for a single emitter, at short time scale (\t\ < 20 ns). Excitation was carried out 
at the maximum available cw power of 9 mW. Integration duration was T = 590 s, 
Ri ~ 37000 counts/s, Rn — 48700 counts/s and timebin w = 0. 17 ns. Left scale: correspond- 
ing normalized intensity correlation function g( 2 )(f). The experimental data are shown as 
red dots, while the sòlid blue line represents a fit by a convolution of Eq.(7) with the meas- 
ured instrumental response function. (b) Evolution of parameter X\ (red dots) as a function 
of the excitation power with linear fitting (in blue), according to Eq.(4). 



The normalization of the g^ correlation function is inferred from the number of coinci- 
dences c(t ) recorded during an integration period T 

^'W-^r-' (9) 

Kl K2 T W 

where the factor R1R2T corresponds to the coincidence rate for an equivalent Poissonian photon 
number distribution, associated with R\ 2 counts rates on each detector of the HBT setup and a 
time bin of width w. Note that this relation between c{t) and g^ 2 \t) neglects any contribution 
from random background light emission, as we can assume from the high value of the signal- 
to-background ratio for the observed emitter shown in Fig. l(a). 

The correlation function g^(t), and the associated measured number of coincidences are 
represented on Fig.3(a) at short time scale. It shows a distinct minimum at zero delay, with a 
residual non-zero value caused by convolution with the instrumental response function (IRF) 
of the detection setup. 

This function is independently measured to be a 1 .2 ns FWHM Gaussian function, by correla- 
tion of Ti:Sa làser 150-fs làser pulses. Taking into account the IRF time, the measured intensity 
correlation function at short time-scale can be well fitted by the convolution of the IRF with the 
g( 2 ^ function given by Eq.(7) for a single emitter. 



By fitting the g^ function curve at short time scale for different excitation powers, we obtain 
the value of X\ dependent on the excitation power (Fig.3(b)). At zero excitation power, X\ is 
equal to the transition rate r2i, which is by definition the invert of the excited-level lifetime. For 
the single emitter studied in Fig.l(a), we get ri\ = 0.47 ns -1 by extrapolating the linear fit to 
zero excitation power. Excited-level lifetime is then calculated to be as short as 2 ns. 

Absorption rate from the ground level is deduced from the excitation intensity dependance 
of Ai, by r\2 = al/hv, where a is the absorption cross-section, / the excitation intensity , and 
V the excitation frequency. From X\ at different excitation intensities, the power-independent 
absorbing cross-section is evaluated to be rj ~ 1.7 x 10~ 16 cm 2 . 




Fig. 4. (a) Number of photon coincidences c(f) (right scale) recorded with the HBT setup 
for a single emitter, on long time-scale (|í| > 20 ns). The excitation power was 9 mW. In- 
tegration duration was T = 605 s, R\ ~ 38600 counts/s, R2 — 33400 counts/s and timebin 
w = 2.3 ns. Left scale: corresponding normalized intensity correlation function g( 2 > (f ). The 
experimental data are shown as red dots, while the sòlid blue line represents a fit by a con- 
volution of Eq.(8) with the measured instrumental response function. The dashed line indi- 
cates the normalization level corresponding to Poissonian photon-number statistics. Note 
that the minimum value of g' 2 ' at t = appears higher than in Fig.3(a) due to choice of a 
larger timebin for histogram plot. (b) Evolution of parameter %i (red dots) as a function of 
the excitation power, with fit (in blue) according to Eq.(5). 



On long time scale, we observe g^ (/) > 1, corresponding to a bunching effect that is due to 
leakage towards the "dark" metastable level. By fitting the g^ (?) function at long time-scale ne- 
glecting the antibunching part, we obtain X2 as a function of the excitation intensity (Fig.4(b)). 
Intersystem crossing rates ^23 and are evaluated by fitting X2 using Eq.(5), together with ru 
and r%\ deduced from the short time-scale measurements. 

Considering r23 as a constant, we observe that the bunching effect depends on the làser 
power, indicating that 687 nm excitation light contributes to deshelving the metastable level. 
Indeed, data are well fitted by taking into account a 1 dependence on the excitation intensity 
(or equivalently on làser excitation power P e x) of the form = r9j (1 + j3P e x), giving ^3 = 
2.75 /is -1 , = 1.71 /is -1 and /3 = 0.102mW _1 . Note that such a dependance of deshelving 
rate with excitation intensity is also observed for single molècules [24]. 

With all parameters obtained above, the detected fluorescence rate is given by 

^ = T7detr? Q7 — -, — -, — — T, (10) 



where R is the sum of count rates on the two APDs in the HBT setup. Parameters Tjdet and ï]q 
represent the overall detection efficiency and the photoluminescence quantum yield, respec- 



tively. With all the decay rates obtained following the fitting procedure of X\ and X2 intensity 
dependence, we obtain ridetlQ = 4.2 x 10~ 4 . 
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Fig. 5. Photoluminescence intensity measurement (plain circles with error bars) corre- 
sponding to the sum of the counting rates on the two detectors of the HBT setup, vs. ex- 
citation làser power. The line is a fit to the data by Eq.(10), taking into account the vàlues 
of rates r mn determined from measurements of parameters X\ and Xi. The fit is therefore 
achieved with a single free parameter left, i.e. the overall efficiency r/det x i?q- 

In our experimental setup, 77^ is limited by the collection efficiency of the microscope objec- 
tive (microscope objective with NA= 0.95 working in air without immersion oil, rj co i ~ 2.6%), 
sphericity aberration (r; ab ~ 20%), transmittance of the objective (Tj tr ans — 80%), transmittance 
of the òptics (Iens, dichroic mirror and füters, íj pt — 35%), and quantum efficiency of the sil- 
icon APDs in the near-infrared region (íJapd — 55%). Then the total detection efficiency is 
TJdet = ^coi x ï7ab x i?trans x TJopt x i?apd — 8 x 10~ 4 . We finally infer a rough estimate of the 
photoluminescence quantum yield as rjg ~ (52 ±20)%. The quantum yield obviously needs 
a more precise evaluation. Our estimation however indicates that the measured short excited- 
level lifetime is presumably associated to strong radiative oscillator strength, and not to efficient 
photoluminescence quenching through fast non-radiative decay processes. 

4. Conclusion 

We have studied, at the single-emitter level, the photophysical properties of a near-infrared 
emitting colour centre, probably of the Nickel-Nitrogen related type. Assuming a three energy 
levels model, we have measured the key parameters associated to its photoluminescence using 
the intensity correlation function both at "short" and "long" time-scale, these two limits be- 
ing compared to the emitter excited-level lifetime. The studied single diamond colour centre 
has very favorable properties for applications in single-photon quantum information processing 
and quantum cryptography. A complementary study will be to determine the optimal excita- 
tion wavelength, in order to design a customized pulsed excitation làser to set up an efficient 
triggered single-photon source. 
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